Estimating the uncertainty in passive-microwave rain-retrievals 
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Introduction 

Current passive-microwave rain-retrieval methods are 
largely based on databases built off-line using cloud models. 
The vertical distribution of hydrometeors within the cloud has 
a large impact on upwelling brightness temperatures ([3], [5]). 
Thus, a forward radiative transfer model can predict off-line 
the radiance associated with different rain scenarios. To esti- 
mate the rain from measured brightness temperatures, one sim- 
ply looks for the rain scenario whose associated radiances are 
closest to the measurements. To understand the uncertainties 
in this process, we first study the dependence of the simulated 
brightness temperatures on different hydrometeor size distri- 
bution (DSD) models. We then analyze the marginal and joint 
distributions of the radiances observed by the Tropical Rain- 
fall Measuring Mission satellite and of those in the databases 
used in the TRMM rain retrievals. We finally calculate the co- 
variances of the rain profiles and brightness temperatures in the 
TRMM passive-microwave database and derive a simple para- 
metric model for the conditional uncertainty, given measured 
radiances. These results are used to characterize the uncertainty 
inherent in the passive-microwave retrieval. 


Effects of the DSD on the radiances 


Most radiative transfer models currently used to calculate 
the expected brightness temperatures ( Tb ) associated with rain 
events assume that the rain drops are distributed according to 
the Marsall-Palmer drop size distribution (DSD) (e.g. [4]), 

N{D)dD = N 0 e~* D dD (1) 


with N 0 = 0.08 cm~ 4 and A = 41 R~ 0/li cm~ l , where R is the 
rain rate. If we calculate the rainfall from this equation, the 
result ( Rpost ) is quite different fron the original rainfall R used 
to determine A. We therefore modified No to make Rpost — 
R. The study has been done with these two Marshall-Palmer 
distributions. A different T-distribution was proposed in [1] : 

N(D) = N 0 {R,S",D") ’ D " ) e - A{ *' S "’ D ”) D , (2) 


with 3 uncorreleted variables: the rain rate /?, the normalized 
mass-weighted mean drop diameter D u and relative deviation 
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Based on the TOGA/COARE data, we found D " = 1.13 ±0.32 
and S" = 0.39 ±0.025. To compare the effets of these distri- 
butions, we first determined 5" and D " in (2) to fit the original 
(MP0) and the adjusted (MP) Marshall-Palmer distribution, by 
minimizing the channels 10.7 GHz and 37.0 GHz. Table I and 
figure I confirm that the Marshall-Palmer distribution with the 
adjusted No is very close to distribution (2). 
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TABLE I 

Parameters which minimize the calculated 10.7 GHz and 37.0 GHz radiances. 


With a homoganaoua atmoaphara 



Fig. 1. Radiative transfer model calculation 

Manifestly, the original and widely used Marshall-Palmer 
distribution implies a larger D", hence larger drops for the same 
brightness temperatures. As figure 1 shows, the 10.7GHz- 
based R-Tb correspondence under-estimates the rainfall by 
about 15% on average, whereas the 37GHz-based correspon- 
dence under-estimates slightly at low rain rates but over- 
estimates subtantially beyond about 8 mm/hr, when scattering 
effects become important. We found that the effect of the DSD 
was minimized near 16 GHz. These results are confirmed by 
cloud model simulations from the TRMM database, as figure 2 
shows. 
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This means that the entire rainfall distribution in the atmo- 
sphere can be described to first order by the vertically averaged 
rain rate R , [ and the constant values of the means of R* ly • * • ,R* 5 . 
When reconstructed in this way, the values for the rain rates 
were within 0.2 % of the original values (see table II). 



rms deviation (mm/hr) 
in mm/hr 

relative deviation 

R 1 

8.95 10- 2 


R 2 

7.89 10“ 2 


R 3 

5.98 10" 2 

1.26953 IQ" 3 

R 4 

4.05 10- 2 


R 5 

1 2.66 10“ z 



TABLE II 

RMS error in the rainrate estimated for each layer from the mean rainrate 
(R*i ) and E {/?2 }*•*•»£ {^5 } 


Fig. 2. Comparison between different DSD calculation. 


Re-parametrization of R 

Assume for simplicity that the typical atmosphere has 5 lay- 
ers, with R constant in each. To understand the joint behavior of 
the rain and the radiances, one must compute their covariance. 
If , /?5 are the 5 eigenvectors of the covariance matrix of 

log(R) calculated for the TRMM cloud simulations database, 



This result, obtained from the TRMM cloud model database, 
is quite similar to the ones found using the TRMM radar data. 
Using five 1-km layers, the eigenvalues for convective events 
in October 1998 were 12.46 > 5.1 > 0.94 > 0.3 > 1. 10~ 2 and 
the orthonormalized eigenvector a, log R l for the first eigen- 
value 12.46 had 0.28 < (a,*) < 0.5, which compares favor- 
ably with our predicted vadue ^ ~ 0.447. So the TRMM 
radar’s R f l is indeed also essentially the average rainfall. For 
stratiform events, the TRMM radar’s eigenvalues were 7.6 > 
1.86 > 0.4 > 0.14 > 8 10 -3 , and the coefficients a { of R f l were 
in the range 0.35 < (a/) <0.51. In both cases, the eigen- 
vector for the second value was (a 1,^2 ><*3,— <* 4 ,— <* 5 )* with 
0.22 < a\ y a 2 ,a 4 ,a$ < 0.52 and a 0.05. Thus, one can 
describe the rainfall to second order using the mean rain rate 
value and a simple difference between the rain in upper layers 
and that in lower layers. 

Estimation of R from Tb 

Since we can describe the vertical rainfall distribution to first 
order using a single variable R f l and constants R f ly • • • y R$, the 
estimation problem is reduced to calculating the conditional 
mean and variance of R f { from the brightness temperatures. To 
find the best relation between R f ] and a combination of the 
TRMM brightness temperatures V at 10.7 GHz, 19.3 GHz, 
21.3 GHz, 37 GHz and 85.5 GHz, we maximized the expec- 
tation E{R / l . T'} keeping E{T' 2 } constant, where: 

r = Vi c r-E{r }) (6) 


we found thatRj ~ -j- Xi log #;• Moreover, the eigenvalues 
are 
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The optimal v, minimize the scatter between T and R 1 . Once 
the v { are found, one can easily compute the mean and variance 
of R* given V . We considered many combinations of passive 
channels, from 5 to 2 brightness temperatures. The results were 
worst when we did not use the 10.7 GHz channel or when we 
used two polarizations of the same frequency. We obtained rea- 
sonably good results using the five vertically polarized chan- 
nels, namely an average R.M.S. uncertainty of 27.0 % on the 
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Fig. 3 . Sample joint distribution, from the TRMM cloud model database, of 
V and the log-average rain rate . 

mean rain rate. Our best results were obtained with horizontal 
polarizations, namely with 

T' = 0.41 7 io. 7 H +0.36 T l93H +0.79 T 2U v 

-0.18 737// —0.182 785.5// (7) 

for which the average R.M.S. uncertainty was 26.4 % and with 

7' = 0.53 7to.7f/ +0.82 7i9.3// -0.2 T 37 // ( 8 ) 

for which the average R.M.S. uncertainty was 26.8 %. 

Figures 3 and 4 illustrate the results obtained with the 7' 
above. Figure 3 shows plots of 7' versus R! v for 7' as in (7) 
and ( 8 ) above and for a sub-optimal 7 ; . Figure 4 shows the 
reconstructed near-surface rain rate R\ plotted against the orig- 
inal. 

Conclusion 

Our study of the joint behavior of the rain in a horizontally 
stratified atmosphere and the associated microwave radiances 
shows that the single most crucial variable characterizing the 
rain profile is the vertically averaged rain rate, followed as 
a distant second by the difference between the high-altitude 
sub-ffeezing-level rain and the precipitation closer to the sur- 
face, the remaining eigen-variables having negligibly small 
variances implying that they can safely be considered constant 
(equal to their respective means). The study also shows that 
a judiciously chosen linear combination of the brightness tem- 
peratures can estimate the rain quite adequately, with an aver- 
age R.M.S. uncertainty (due to the variations accounted for in 
the TRMM cloud model database) of about 27%. The DSD 
does affect the brightness temperatures, and hence the eventual 
retrievals. Below 16 GHz, where scattering is not significant, 
the lager the mean drop size the smaller the rain associated with 



Fig. 4 . (found from V) versus the original using 3 horizontal brightness 

temperatures (7 io. 7 ,79.3 and 737). 

a given brightness temperature. The effect is reversed above 16 

GHz for higher rain rates. 
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